Effect of aperture field distribution on the maximum radiated power at atmospheric pressure
Zhao Pengcheng, Guo Lixin
School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China

 

† Corresponding author. E-mail: pczhao@xidian.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 61501358, 11622542, 61431010, and 61627901) and the Fundamental Research Funds for the Central Universities, China.

Abstract

The air breakdown in the high-power antenna near-field region limits the enhancement of the radiated power. A model coupling the field equivalent principle and the electron number density equation is presented to study the breakdown process in the near-field region of the circular aperture antenna at atmospheric pressure. Simulation results show that, although the electric field in the near-field region is nonuniform, the electron diffusion has small influence on the breakdown process when the initial electron number density is uniform in space. The field magnitude distribution on the aperture plays an important role in the maximum radiated power above which the air breakdown occurs. The maximum radiated power also depends on the phase difference of the fields at the center and edge of the aperture, especially for the uniform field magnitude distribution.

1. Introduction

The high-power microwave has important applications in many fields, such as the long-distance wireless energy transport, creation of ionospheric and atmospheric ducts for long-range communications, high-power radar systems, and plasma chemistry.[1] In recent years, the electric field radiated from the aperture antenna reaches several MV/m that can exceed the breakdown electric field at atmospheric pressure.[2,3] Once the air breakdown occurs, the attendant plasma strongly hinders the microwave radiation. This suggests that the radiated power reaches the maximum (or critical) value when the maximum electric field in the near-field region is equal to the breakdown electric field. Therefore, it is very important to study and understand the maximum radiated power limited by the air breakdown.

The air breakdown caused by the high-power microwave has been investigated by many scholars in recent years.[49] Hidaka et al. observed the variation of the plasma pattern produced in the air breakdown with the pressure, which was well reproduced and explained by the theory model of Boeuf et al.[4,5] Experiments of Yang et al. showed that the breakdown threshold of the short-pulse microwave depends strongly on the gas pressure and species.[6] The electron fluid model with the accurate rate coefficients was employed by Zhao et al. to simulate and reveal the air breakdown caused by the high-power microwave.[7,8] In general, we do not desire the air breakdown to occur, and therefore determining the maximum radiated power above which the air breakdown occurs is important to avoid the breakdown. The field equivalent principle and the empirical formula of the ionization rate were employed by Zhang et al. to study the maximum radiated power of the aperture antenna, where the field phase on the aperture was assumed to be uniform.[9] However, the field phase on the aperture can be nonuniform, which has an important impact on the pattern of the radiated field. It can also be expected that the maximum radiated power changes with the field phase distribution.

In this paper, we study the air breakdown process in the near-field region of the circular aperture antenna at atmospheric pressure using a model coupling the field equivalent principle with the electron number density equation. The electron diffusion that may have an influence on the breakdown process in the nonuniform field is included in the electron number density equation. We first simulate the evolution of the electron number density in the breakdown process, and show how the electron ionization, attachment, and diffusion affect the evolution. Then we focus on the dependence of the maximum radiated power on the field phase and magnitude distribution on the aperture.

2. Model

The plasma produced in the air breakdown process has little influence on the incident microwave until its number density reaches the critical value , where ε0 is the permittivity of free space, ω is the angular frequency of the microwave, and me and qe denote the charge and mass of electron, respectively.[10] In this paper, we focus only on the process before air breakdown and predict the maximum power radiated from the circular aperture antenna above which the air breakdown occurs. We first adopt the field equivalent principle to solve the radiated electric field in the near-field region, which can be approximated as a free space. Then we use the electron number density equation to describe the evolution of the electron number density in the breakdown process. Finally, the definition of the maximum radiated power is presented.

2.1. Field equivalence principle

The equivalent electric current on the aperture alone can be written as[11] where Ea is the electric field on the aperture and η0 is the wave impedance of the free space. The vector potential corresponding to the electric current is given by where is the wave number, μ0 is the permeability of the free space, and R is the distance between the field point r and the source point r′. By using Eq. (2), the radiated fields can be obtained as follows:

2.2. Electron number density equation

Considering the fact that the radiated electric field in the near-region is nonuniform, the electron diffusion term is included in the electron number density equation. Since we focus on the process in which the electron number density grows from the initial level to the critical value above which the incident wave is disturbed, the electron number density is low and then the electron–ion recombination can be ignored. Therefore, the electron number density can be expressed as[12] where ne denotes the electron number density, νi and νa are the ionization rate and attachment rate, respectively, and Deff is the effective electron diffusion coefficient.

The νi and νa are obtained from the following empirical scaling laws:[12] where p is the air pressure, Ec ≈ 3200p is the critical field at which the ionization balances the attachment, and Eeff is the effective electric field. Note that equation (6) holds only when Eeff/Ec lies between 1 and 3. Eeff can be written as where the collision frequency νm is approximated as Deff is expressed as[5] where α is the local Maxwell relaxation time, μe = qe/meνm and μi = μe/200 are the electron mobility and ion mobility, respectively, Da = (μi/μe) De is the ambipolar diffusion coefficient, De = μeTe/e is the free diffusion coefficient, and the electron temperature Te is approximated as[12] Note that equation (12) is valid only when Eeff/p ≤ 100 V⋅cm⋅Torr−1.

2.3. Maximum radiated power

When the electron diffusion term can be ignored, integration of Eq. (5) yields where ne0 is the initial electron number density. The electron number density grows exponentially and the breakdown is defined to occur when the density has increased to the critical level above which the microwave is disturbed obviously by the breakdown plasma.[5] Under the breakdown electric field, equation (13) can be rewritten as where tp is the microwave pulse width. By introducing Eqs. (6)–(8) into Eq. (14), the rms breakdown electric field is obtained as follows: Let tp → ∞ in Eq. (15), the rms breakdown electric field for the continuous microwave can expressed as When the effect of the electron diffusion on the breakdown process is obvious, we must adopt Eq. (5) to determine Eb.

The power density corresponding to the breakdown electric field is where ‘*’ represent the conjugate operation, and |E| = Eb. When the maximum radiated power density reaches Pb, the antenna radiated power approaches the upper limit where γmax is the ratio of the maximum radiated electric field to the maximum aperture electric field, is the normalized electric field on the aperture, and S′ is the aperture. It is clear that γmax ≥ 1.

3. Results and discussion

The air pressure is taken as 760 Torr that corresponds to one atmospheric pressure. At the initial time, the electron number density ne0 = 1 × 106 m−3 is uniform in the near-field region. The circular aperture with the radius of 0.2 m is placed on the xy plane, whose center is located at the origin of the coordinate system. It is assumed that the microwave electric field on the aperture has one component, i.e., Ea = exExa. When the field distribution on the aperture is taken to be cylindrical symmetric, the field radiated from the aperture is also cylindrical symmetric. Therefore, the three-dimensional problem can be reduced into a two-dimensional one and we can consider only the electric field and the electron number density on the yoz plane. We focus on the breakdown prediction of the continuous microwave with a frequency of 2.85 GHz. The microwave is a monochromatic wave, whose magnetic and electric field can be directly solved by Eqs. (3) and (4).

Figure 1 shows the electric field radiated from the aperture on which the phase and magnitude of the field are uniform. It is obvious that the radiated field is nonuniform. When the maximum radiated field is taken to be slightly higher than the critical field Ec, for example, 1.1Ec, the evolution of the electron number density is simulated, as shown in Figs. 2 and 3. In the high field region, the ionization rate is larger than the rate of electron loss due to the attachment and diffusion, so the electron number density at 50 ns in Fig. 2(a) shows an obvious increase compared with the initial one. We also find that, with further increase in time, the electron number density in the high field region becomes higher, while its profile is similar to the one at t = 50 ns, as shown in Fig. 2(b). This is because the electron diffusion is much less than the electron ionization and attachment, as shown in Fig. 3. The phenomenon is also observed at other aperture field distributions. The secondary electron emission is not considered. This is because at a high pressure the mean secondary emission yield is much lower than unity, and the ionization becomes the dominant mechanism to generate electrons.[13]

Fig. 1. (color online) The ratio of radiated electric field to maximum aperture electric field. The phase and magnitude of the field on the aperture are uniform.
Fig. 2. (color online) The spatial distribution of the electron number density at t = 50 ns (a) and 200 ns (b) under the radiated field corresponding to Fig. 1. The maximum radiated electric field is equal to 1.1Ec.
Fig. 3. (color online) Comparison among electron ionization, attachment, and diffusion at the point where the radiated field is maximum. The radiated field profile corresponds to Fig. 1 and its maximum value is equal to 1.1Ec.

We consider three different distributions of the aperture field magnitude, which are the uniform, parabolic taper with pedestal, and parabolic taper distributions, respectively.[11] The distributions of parabolic taper with pedestal and parabolic taper can be written as follows: The field phase on the aperture is assumed as where ϕmax is the phase difference of the fields at the center and edge of the aperture. Figure 4 shows the ratio γmax of the maximum radiated electric field to the maximum aperture electric field for the three different aperture field distributions. It can be seen from this figure that the field magnitude distribution plays an important role on γmax. Of the three field magnitude distributions, γmax corresponding to the uniform field magnitude distribution is highest. γmax for the uniform field magnitude distribution depends obviously on the phase difference ϕmax, while for the other two field magnitude distributions, the dependence appears only at ϕmax < π.

Fig. 4. (color online) The ratio of maximum radiated electric field to maximum aperture electric field for three different aperture field distributions.

Figure 5 shows the maximum radiated power above which the air breakdown occurs. The aperture field distributions can refer to Eqs. (19)–(21). Since the electron diffusion has small influence on the breakdown process (see Fig. 3), equation (16) can be adopted to predict the breakdown electric field Eb. Under the fixed air pressure and field magnitude distribution, Eb remains unchanged, and (see Eq. (18)). Therefore, the variation of Pmax with ϕmax has a contrary trend with that of γmax, as shown in Figs. 4 and 5. Of the three different field magnitude distributions, both γmax and corresponding to the uniform field magnitude distribution are highest. As a result, the corresponding Pmax is highest around ϕmax = 2.5π.

Fig. 5. (color online) The maximum radiated powers for three different aperture field distributions.

To validate the breakdown prediction, the breakdown thresholds obtained from Eq. (15) are compared with the experiment results of Cook et al. (see Ref. [14]) under f = 110 GHz and tp = 2500 ns, as shown in Fig. 6. It can be seen from this figure that the breakdown prediction agrees very well the experimental data. The fluctuation in the experimental data is caused by the seed electrons.[2]

Fig. 6. (color online) Breakdown thresholds obtained from Eq. (15) and experiments of Cook et al. under f = 110 GHz and tp = 2500 ns.
4. Conclusion

The model coupling the field equivalent principle and the electron number density equation is presented to study the air breakdown process in the near-field region of the circular aperture antenna, and to predict the maximum radiated power above which the air breakdown occurs. The terms of the electron ionization, attachment, and diffusion are included in the electron number density equation. The magnitude of the aperture field is assumed to have the uniform, parabolic taper with pedestal, and parabolic taper field distributions, respectively, and its phase is nonuniform. The results show that, although the electric field in the near-field region is nonuniform, the electron diffusion has small influence on the breakdown process when the initial electron number density is uniform in space. The distribution of the aperture field magnitude plays an important role in the maximum radiated power. The maximum radiated power for the uniform field magnitude distribution depends on the phase difference of the fields at the center and edge of the aperture, while for the other two field magnitude distributions, the dependence appears only at small phase differences.

Reference
[1] Golubev A I Sysoeva T G Terekhin V A Tikhonchuk V T Altgilbers L L 2000 IEEE Trans. Plasma Sci. 28 303
[2] Ford P J Beeson S R Krompholz H G Neuber A A 2012 Phys. Plasmas 19 073503
[3] Chang C Verboncoeur J Wei F L Xie J L Sun J Liu Y S Liu C L Wu C 2017 IEEE Trans. Dielectr. Electr. Insul. 24 375
[4] Hidaka Y Choi E M Mastovsky I Shapiro M A Sirigiri J R Temkin R J 2008 Phys. Rev. Lett. 100 035003
[5] Boeuf J P Chaudhury B Zhu G Q 2010 Phys. Rev. Lett. 104 015002
[6] Yang Y Yuan C Qian B 2012 Phys. Plasmas 19 122101
[7] Zhao P Guo L Li H 2015 Chin. Phys. 24 105102
[8] Zhao P Guo L Shu P 2017 Chin. Phys. 26 029201
[9] Zhang J Wang J 2011 IEEE Electromagn. Compat. 53 540
[10] Zhao P Liao C Feng J 2015 Chin. Phys. 24 025101
[11] Stutzman W L Thiele G A 1998 Antenna Theory and Design 2 New York Wiley 275 283
[12] Zhu G Q Boeuf J P Chaudhury B 2011 Plasma Sources Sci. Technol. 20 035007
[13] Kim H C Verboncoeur J P 2006 Phys. Plasmas 13 123506
[14] Cook A Shapiro M Temkin R 2010 Appl. Phys. Lett. 97 011504